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Abstract 

In this Letter we present first principle calculation of surface energies of rock-salt (Bl), zinc- 
blende (B3), and wurtzite (B4) A1N allotropes. Out of several low- index facets, the highest 
energies are obtained for mono-atomic surfaces (i.e. only by either Al or N atoms): lf\n\ — 

410 me V/ A 2 , 7^ = 346meV/A 2 , 7^ = 360meV/A 2 , and 7 { B 4 001} = 365meV/A 2 . The 
difference between Al- and N-terminated surfaces in these cases is less then 20meV/A . 
The stoichiometric facets have energies lower by lOOmeV/A or more. The obtained trends 
could be rationalised by a simple nearest-neighbour broken-bond model. 
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1. Introduction 

Aluminium nitride is a material with a fascinating variety of applications spanning from 
optoelectronic [l| and acoustic [2] devices to the beneficial influence on mechanical proper- 
ties (e.g., increased hardness) and performance of protective hard coatings It crystallises 
in the wurtzite structure (B4, space group PQ^mc). Under non-equilibrium deposition pa- 
rameters, its metastable zinc-blende variant^ (B3, space group F43m) can be stabilised j4J. 
Finally, under pressures above « 16GPa it transforms to the rock-salt structure (Bl, 
space group Fm3m). 

The knowledge of surface energies and the corresponding trends is important in several 
areas. As discussed by Gall et al. [6] for the case of TiN, surface energy is a decisive 
parameter for the thin film microstructure when grown under near-to-equilibrium conditions, 
i.e. when thermodynamics rather than kinetics control the texture formation. This is the 
case, for example, of the chemical vapour deposition of A1N for optoelectronic applications. 
Equally important is this quantity for the discussion of decomposition pathways of, e.g., 
unstable Ti^^AL-N 0, 0]. Another emerging area of interest is the stabilisation of rock- 
salt A1N using a multilayer architecture for Ti 1 _ a .Al, r N/AlN or Crj^Al^N/AlN coatings j^j, 
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as it strengthens the material, and has the potential for stopping and/or deflecting cracks 
flpj ] . Due to the nanostructured character of these composites, surface energy represents a 
non-negligible contribution for the overall energetic balance 11 



Although there exist plenty of papers, both on calculations as well as on experimental 
results, a systematic study of the surface energies is missing. In this Letter we use quantum 
mechanical calculations to obtain energies of several low-index clean surfaces (i.e., no surface 
reconstruction is considered here) of the Bl, B3, and B4 allotropes. Since we cover several 
surface orientations as well as crystallographic structures, these results represent a coherent 
data set which can be directly used as is in above described applications. 



2. Calculational details 

Quantum mechanical calculations within the framework of Density Functional Theory 
(DFT) were performed using the Vienna Ab initio Simulation Package (VASP) [12||.The pro- 
jector augmented wave method pseudopotentials used Generalised Gradient Approximation 



(GGA) [13j as parametrised by Perdew and Wang 14]. The first Brillouin zone (1BZ) of the 
bulk material unit cells was meshed with 5x5x5, 8x8x8, and 9x9x6 Appoints distributed 
according to the Monkhorst-Pack scheme for the Bl, B3, and B4 allotrope, respectively. The 
sampling was adopted for the surface supercells (see below), with only 1 fc-point along the 
direction perpendicular to the surface. The total energy was obtained by integration over 
the whole 1BZ using the tetrahedron method. These parameters together with plane wave 
cut-off energy of 450 eV ensure accuracy of the total energy of ~ lmeV/atom. 

For each structure and each surface orientation, an iV ce ii-atom primitive cell was con- 
structed such that the particular facet was perpendicular to the z-axis. Subsequently, n 
primitive cells were stacked along the z-axis followed /i vac thick vacuum. Periodic boundary 
conditions were applied on those slabs. If the surface area is A and contains iV sur f atoms, 
the surface energy is calculated as 

7(n, h vac ) = ~-r(E 8 l ah - E hu i k (nN ceU + iV sm . f )) (1) 

where -Esiab is the total energy of the slab while -Ebuik is the bulk total energy per atom of 
the respective structure. Positions of all atoms were relaxed during calculations of -E s i a b- 

Finally, 7(71, h vac ) is converged simultaneously with respect to the number of primitive 
cells in the slab, n, and the vacuum thickness, h vac , in order to rule out undesired surface 
interactions either though the bulk material or though the separating vacuum. 



3. Results and discussion 

The calculated bulk lattice parameters (see Tab. [1]) agree well with both calculated and 
experimental values previously published. The lowest energy of formation is obtained for the 
wurztite B4 structure, closely followed by the zinc-blende B3 (AEf 3 ~ m = 0.0213 eV/at.). 
The rock-salt structure Bl possesses highest energy of formation (AE® 1 ' 34 = 0.1813eV/at.). 
These values correspond well with data calculated by Siegel et al. |l5| . 
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Table 1: Lattice constants, a and c, total energy, -Etotj and energy of formation, Ef, of individual A1N 
structures. Calculated and experimental lattice parameters from the literature are given for comparison. 
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Figure 1: Convergence of the Al-terminated {111} surface energy of B1-A1N. 



In this study we consider 4 low-index surfaces ({100}, {110}, {111}, and {112}) for the 
cubic Bl and B3 structures, and 3 different orientations ({0001}, { 1 100}, and {1120}) for 
the hexagonal B4 phase. The respective primitive cells contain 4, 4, 6, and 12 atoms for 
cubic structures and 4, 8, and 8 for the wurztite one. Their exact description is given in 



Appendix A\ It is worth noting that the Bl-{111}, B3-{100}, B3-{111}, and B4-{0001} 



planes contain only atoms of one specie. Consequently, in these cases we have calculated 
separately Al- and N-terminated surfaces. 

A typical example of the surface energy as a function of the slab and vacuum thickness is 
shown in Fig. [TJfor the Al-terminated Bl-{111} surface. As in all other cases, the difference 
between results for 10 A and 15 A is negligible, while 5 A gives different 7 values. The slab 
thickness of ~ 20 A is sufficient in this case, but for other orientations thicker slabs are 
needed. 

The resulting surface energies are shown in Fig. [2J The highest values of 7 (between 
340-420 meV/A are for all structures obtained for mono-specie surfaces. Additionally, the 
difference between surfaces terminated by Al or N atom are ~ 20 meV or less. This means 
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that in the cases when only the orientation rather than the exact structure of the surface 
is known, one can used the average value of Al- and N-terminated surfaces as a reasonable 
estimate without changing e.g., the energetic order of individual facets. The surface energy 
of facets formed by both atomic species (stoichiometric facets) is by at least lOOmeV/A or 
more lower than the mono-atomic ones. The lowest surface energies of individual structures 
posses {100} planes for Bl, {110} planes for B3 and {1100} planes for B4. 

The observed trends can be rationalised by a simple nearest-neighbour broken-bond 



model [18|, [19j. If one assumes that the main contribution to the surface energy is the 
"penalty" related to breaking bonds that cross the actual surface, then 7 should be pro- 
portional to AbondsM, where Ab on ds is the number of broken bonds per surface area A of 
the primitive cell. This density of broken bonds is also plotted in Fig. |2j It captures the 
surface energy trends surprisingly well, the only deficiencies being the mutual relations of 
the cubic {110} and {112} facets. Since no plane can cut more then 3 bonds of the oc- 
tahedrally coordinated sites in the Bl structure (which is the case for {112} plane), and 
since the surface area per atom increases for planes with higher Miller indices, the nearest- 
neighbour broken-bond model suggests that the surface energy of all those planes is less 
than 7{n2} = 270meV/A . The situation is slightly more complex for the octahedrally 
coordinates B3 and B4 structures. If always a pair of mirror interface is considered (as in 
the case of this study where e.g., 7^0001} = \ (ifmoi) + 7^iooi)))' then the maximum number 
of broken bonds per surface atom is 2. The only possibility how to get a surface with atoms 
having only 1 broken bond is that 2 out of 4 bonds lie in the surface. In all other cases, the 
surface atoms have 2 broken bonds per atom on average. Out of those, the B3-{111} and 
B4-{0001} planes have the smallest area per atom, thus justifying why 7^i 3 n} and 7{ 3 ooi} are 
the highest surface energies for those structures. 

4. Conclusions 

We used the quantum mechanical calculations to obtain surface energies of low-index 
surface facets of rock-salt Bl, zinc-blende B3, and wurtzite B4 allotropes of A1N. No surface 
reconstruction was taken into account. The highest surface energies were obtained for mono- 
specie planes i.e., 7{m} = 410 me V/A' , 7fioo} = 346 me V/A , 7{m} = 360meV/A , and 

^foooi} = 365meV/A . The observed trends could be rationalised by the nearest-neighbour 
broken-bond model: the higher the areal density of broken bonds, the higher the surface 
energy. 

Acknowledgements 

Financial support by Christian Doppler Research Association as well as by the START 
Program (Y371) of the Austrian Science Fund (FWF) is greatly acknowledged. 

References 



[1] S. C. Jain, M. Willander, J. Narayan, R. V. Overstraeten, J Appl Phys 87 (2000) 965. 

4 



Appendix A. Structures 



Table A. 2: Lattice vectors and fractional coordinates of atoms for the primitive cells used for construction 
of {100}, {110}, {HI}, and {112} surfaces of the Bl structure. The N-terminated {111} slab is obtained 
by exchanging Al and N lattice sites. 
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Table A. 3: Lattice vectors and fractional coordinates of atoms for the primitive cells used for construction 
of {100}, {110}, {111}, and {112} surfaces of the B3 structure. The N-terminated {100} and {111} slabs 
are obtained by exchanging Al and N lattice sites. 
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Figure 2: Surface energies for low-index planes of A1N allotropes. Red curves represent calculated density 
of broken bonds as discussed in the text. 



Table A. 4: Lattice vectors and fractional coordinates of atoms for the primitive cells used for construction 
of {0001}, {1100}, and {1120} surfaces of the B4 structure. The N-terminated {0001} slab is obtained by 
exchanging Al and N lattice sites, u is the internal parameter describing shift of N atoms above Al sites 
(wain = 0.382). c and a are the lattice parameters. 





{0001} 


{1100} 


{1120} 


ai 


1/2 


-VS/2 





1 








3/2 








a 2 


1/2 


V3/2 








c/a 








c/a 





a 3 








c/a 








Vs 








V3/2 


Al 


1/3 


2/3 























Al 


2/3 


1/3 


1/2 


1/2 


1/2 


1/6 


1/3 


1/2 





Al 








1/2 





1/2 


1/2 





1/2 


Al 











1/2 


2/3 


5/6 


1/2 


1/2 


N 


1/3 


2/3 


u 





u 








;/ 





N 


2/3 


1/3 1/2 


+ u 


1/2 


1/2 + u 


1/6 


1/3 


1/2 + u 





N 








1/2 


u 


1/2 


1/2 


u 


1/2 


N 











1/2 + u 


2/3 


5/6 


1/2 + u 


1/2 



[2] R. Turner, P. Fuierer, R. Newnham, T. Shrout, Appl Acoust 41 (1994) 299-324. 

[3] P. H. Mayrhofer, A. Horling, L. Karlsson, J. Sjolen, T. Larsson, C. Mitterer, L. Hultman, Appl Phys 
Lett 83 (2003) 2049-2051. 

[4] I. Petrov, E. Mojab, R. C. Powell, J. E. Greene, L. Hultman, J.-E. Sundgren, Appl Phys Lett 60 (1992) 
2491. 

[5] D. Holec, F. Rovere, P. H. Mayrhofer, P. B. Barna, Scr Mater 62 (2010) 349-352. 

[6] D. Gall, S. Kodambaka, M. A. Wall, I. Petrov, J. E. Greene, J Appl Phys 93 (2003) 9086. 

[7] P. Mayrhofer, F. Fischer, H. Bohm, C. Mitterer, J. Schneider, Acta Mater 55 (2007) 1441-1446. 



6 



[8] R. Rachbauer, S. Massl, E. Stcrgar, D. Holec, D. Kiener, J. Keckes, J. Patscheider, M. Sticfel, H. Leitner, 
P. H. Mayrhofer, J Appl Phys 110 (2011) 023515. 

[9] V. Chawla, D. Holec, P. H. Mayrhofer, Comput Mater Sci 55 (2012) 211-216. 

[10] M. Schlogl, J. Paulitsch, J. Keckes, P. H. Mayrhofer, Surf Coat Tcchnol under revi (2012). 

[11] V. Chawla, D. Holec, P. H. Mayrhofer (2012). 

[12] G. Kresse, J. Furthmiiller, Phys Rev B 54 (1996) 11169-11186. 

[13] G. Kresse, D. Joubert, Phys Rev B 59 (1999) 1758-1775. 

[14] Y. Wang, J. P. Perdew, Phys Rev B 44 (1991) 13298-13307. 

[15] A. Siegel, K. Parlinski, U. Wdowik, Phys Rev B 74 (2006) 104116. 

[16] Powder Diffraction File 00-046-1200, International Center for Diffraction Data, PDF-2/Release 2007, 
2007. 

[17] O. Madclung, Semiconductors: Data Handbook, volume 8, Springer- Verlag Berlin Heidelberg, 2004. 

[18] E. Kozeschnik, I. Holzer, B. Sonderegger, J Phase Equilib Diffus 28 (2007) 64-71. 

[19] R. Becker, Ann Phys 424 (1938) 128-140. 



7 



